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Abstract 

We study the coupling of massive fermions to the quantum mechanical dynamics of space- 
time emerging from the spinfoam approach in three dimensions. We first recall the classical 
theory before constructing a spinfoam model of quantum gravity coupled to spinors. The 
technique used is based on a finite expansion in inverse fermion masses leading to the com- 
putation of the vacuum to vacuum transition amplitude of the theory. The path integral 
is derived as a sum over closed fermionic loops wrapping around the spinfoam. The effects 
of quantum torsion are realised as a modification of the intertwining operators assigned to 
the edges of the two-complex, in accordance with loop quantum gravity. The creation of 
non-trivial curvature is modelled by a modification of the pure gravity vertex amplitudes. 
The appendix contains a review of the geometrical and algebraic structures underlying the 
classical coupling of fermions to three dimensional gravity. 

CPT-P71-2006 

1 Introduction 

The subtle interplay between matter and the classical geometry of spacetime is probably the 
deepest physical concept emerging from Einstein's equations. The becoming of this relationship 
in the Planckian regime, where the quantum mechanical aspects of spacetime become predominant, 
is a core issue in any attempt to the quantisation of gravity. In this paper, we address the problem 
of fermionic couplings to three dimensional quantum gravity in the spinfoam approach. 

Spinfoam models [T] are discretized path integrals concretely implementing the Misner-Hawking 
three-geometry to three-geometry transition amplitude [2]. They can be realised from different 
approaches. First, a spinfoam can be seen as a spacetime history interpolating between two canon- 
ical states, i.e as a representation of the physical scalar product [3] of loop quantum gravity [4j. 
More generally, spinfoam models can also appear as a background independent, non-perturbative 
covariant quantisation of BF-like field theories by prescribing a regularization of the path integral 
of the theories combining (random) lattice field theory and group representation theory techniques 
in various dimensions. If the dimensionality of spacetime is different from three, where gravity 
is a BP theory, the local degrees of freedom of general relativity can by inserted directly into 
the discretized path integral by constraining BP theory down to gravity following Plebanski's 
prescription [5] directly at the quantum level [5], [7]. The path integral for gravity, called the 
Barrett-Crane model, is therefore obtained as a sum over a restricted subset of configurations. To 
obtain fully background independent models and to restore the infinite number of degrees of free- 
dom of the theories, one works with a dual field theory formulation [5] generalising matrix models 
to higher dimensions. The idea is to write a field theory over a group manifold (GPT) whose 

* fairbairn@cpt.univ-mrs.fr 

t Unite mixte de recherche (UMR 6207) du CNRS et des Universites Aix-Marseille I, Aix-Marseille II et du Sud 
Toulon- Var. Laboratoire affilie a la FRUNAM (FR 2291). 



1 



Feynman amplitudes yield the spinfoam partition function defined on the cellular decompositions 
dual to the Feynman diagrams of the field theory. This approach has recently been rethought as 
a third quantised version of gravity, including a sum over topologies [HI- 

In this paper, we explicitly construct a spinfoam model describing the coupling of massive 
fermionic fields to three dimensional Riemannian quantum gravity. We derive the spinfoam model 
from a classical action principle by giving a precise meaning to the vacuum to vacuum transition 
amplitude of the theory 



by regularising, in particular, the functional measures. We proceed by expanding the fermionic 
determinant, obtained by integrating out the fermionic degrees of freedom, into closed fermionic 
loops. 

The study of matter couplings to gravitation in the spinfoam approach has recently raised 
a large number of proposals (for models from the canonical perspective see |10j). Two classes 
of models seem to emerge from the overall picture: ideas converging towards unification, where 
matter is somehow already hidden in the models, and constructions where matter is added by 
coupling to the gravitational field. In [TT], a unified picture has been suggested where matter 
emerges as the low energy limit of the models derived from the configurations left out of the path 
integral for gravity when constraining BF theory down to gravity. The same author |12| proposes 
an alternative picture where matter appears as the conical singularities of the triangulated pseudo- 
manifolds arising in the GFT Feynman expansion. In [13j . the matter degrees of freedom are 
encoded in a subset of the finite dimensional irreducible representations of the symmetry group 
of gravity. A supersymmetric three dimensional model based on the representation theory of 
OSp{l I 2) has been derived in [H]. In this approach, fermionic matter is "hidden" in the gauge 
field hence offering another unified picture. From the coupling perspective, different models for 
matter and different techniques have been implemented. Yang-Mills fields have been coupled to 
four dimensional quantum gravity in [15j and [16| using lattice field theory techniques where the 
lattice geometry is determined by the Barrett-Crane model. A spinfoam model for pure Yang-Mills 
has been derived in [T7]. Background independent point particle couplings have been introduced 
in three dimensional quantum gravity [IH] (see [12] for the canonical aspect) by using gauge 
fixing techniques or, equivalently, the de-Sousa Gerbert [20] algebraic description of the degrees 
of freedom of a point particle. GFT models have then followed [H] defining matter couplings to 
a third quantised version of gravity. Recently, a subtle relationship between quantum field theory 
Feynman diagrams and three dimensional spinfoam models has been derived. Feynman diagrams 
have been shown [52], [53j to yield natural observables of spinfoam models. Also, in the no- 
gravity limit of spinfoam models coupled to point particles [24j , one recovers the Feynman diagram 
amplitudes of quantum field theory. A generalisation of the above results to four dimensions has 
recently been proposed [5S] from the canonical perspective. It is shown how the natural higher 
dimensional extension of three dimensional quantum gravity coupled to point particles leads to 
quantum BF theory coupled to strings and branes. The derivation of the dynamics of this proposal 
is currently under study [1^ . Finally, it has recently been realised [27] that the coupling of Wilson 
lines to the de-Sitter gauge theory formulation of gravity leads to the description of the dynamics 
of particles coupled to the gravitational field. 

Here, we concentrate on the interplay between spinor fields and three dimensional quantum 
geometry. The paper is organised as follows. Section two is a review of the classical coupling 
of fermions to the gravitational field. We recall, in particular, the effects of fermionic sources 
on spacetime torsion and curvature. The third section describes the construction of the model. 
We define the manifold and field discretizations, provide a simplicial action, define the fermionic 
loop expansion, and finally perform the integrations in the discretized path integral leading to a 
precise definition of the model. The appendix describes the geometrical and algebraic framework 
underlying the coupling of spinors to classical three-dimensional gravity. 
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2 Classical theory 



In this first section, wc derive the dynamical effects of fermions on the geometry of spacetime. 
All necessary definitions and conventions are given in the appendix. We first introduce the action 
governing the dynamics of the gravitational field before adding the fcrmionic term. We then derive 
the equations of motion of the coupled system. 

2.1 Gravitational and fermionic action principles 

Consider a connected, oriented, compact, three dimensional differential manifold AI endowed with 
an Euclidean metric g with diagonal form rj as our spacetime. We will denote V the bundle of 
g-orthonormal frames, that is, the principal bundle with base manifold M and structure group 
G ~ SO{rj) ~ 5*0(3). Noting V — M'^, we call (tt, the vectorial (fundamental) representation 
of G. Let G = Spin{i]) = Spin{3) denote the spin group associated to the Euclidean metric rj. 
We will choose as a basis of the real Lie algebra g = spin(3) (resp. g = so(3)) the set of generators 
{-^a}a (resp. {Ta}a), CL = 1,2,3. We will assume that the manifold M is endowed with a given 
spin structure, i.e. a G-principal bundle V mapped with a two-to-one bundle homomorphism x 
onto the bundle of orthonormal frames V. Note that the principal bundle V is necessarily trivial. 
Let Lo denote a metric (but not necessarily Riemannian) connection on V image under the Lie 
algebra isomorphism = tt, : 5pin(3) so(3); Xa Ta, of a spin connection tj on V. The data 
of the manifold M together with the non-Riemannian metric connection cu defines the Riemann- 
Cartan structure underlying the framework necessary to model the interaction of fermions with 
the dynamical spacetime M. 

We will use the fact that in the adjoint representation of the Lie algebra spin(3) the matrix 
elements of the images of the generators are given by ^T^,{Xa)^ j = j ~ j, where the symbols 
el J denote the structure constants of the Lie algebra normahsed such that £123 = —£132 = 1- We 
will furthermore use the isomorphism of vector spaces between sptn(3) and y = R"^, regarded as 
the adjoint representation space of 5pin(3), to make no distinction between the Lie algebra indices 
a, 6 = 1,2,3 and the vector space indexes /, J = 1, 2, 3. 

The basic dynamical fields of Einstein-Cartan gravity are the soldering one-form e ~ ® ej 
and the pull-back to M by local trivialising sections of a principal metric connection on V. In 
a local chart {U C Af, : U — ^ M"^), the co-frame decomposes into a coordinate basis of the 
co-tangent space eJ = e^da;^, and the metric connection on M is given by lu^'^ = uij/dx'^ = 
w° (g) T^'^dx^ = -t^jiUJ^dx^. 

Since we are focusing on the inclusion of fermionic matter, it will be convenient to work in 
the spinor representation of 2 -I- 1 gravity (see e.g. [13 )• To this aim, we define the linear map 
a : V ~> End (C^), mapping the Clifford algebra C(3,0) = C(3) (in which the (dual of the) 
vector space V is embedded) onto the Pauli algebra, i.e. the endomorphism algebra of the two- 
dimensional complex vector space V = (see appendix). We consider as the dynamical fields 
of the theory the image e of the soldering form e under the map a and (the pull-back to M 
by a global trivialising section of) a principal connection on the spin bimdlc V, that is, a spin 
connection uj: 

1 i 

e := die) = eJaj and lu = ——cu^'^aiaj ~ —LU^ai, (2) 

where the symbols cr/ = cr(ei) denote the Pauli matrices. 

Hence, we can think of the soldering form e : TpM Af2(C), forall p in M, as a local 
isomorphism mapping any tangent vector into a two by two complex matrix, image in spin space 
of the associated vector in inertial space V: Vw S TpM, e{v) = v = v'aj. More precisely, the 
associated matrix v is traceless and Hermitian. Accordingly, we can think of the co-frame as 
taking value in i5pin(3) ~ isu(2). 
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The torsion T = a{T^ei) and curvature R ~ —jR^'^ajaj of the spin connection are then 
extracted from to torsion and curvature of the metric connection to: 

T' = de^ + e^w-^Ae^ (3) 
R' = du' + A w-^, 

where to' = —^e^ jj^uj''^. 

Gravitational action. The classical dynamics of the three dimensional spacetime M can be 
expressed in terms of the Palatini action evaluated in the spinor representation of spin(3) 

SGR[e,u;]^- ( Tr(eAR[Cj]), (4) 

where Tr is the trace on A'/2(C), and k is related to the Newton constant G by k = SttG. 

Equivalently, we can write the Palatini action in the vector representation of the spin group 

SgrVM f ^UKc' A i?^^M = -- / e/ ARi[uj], (5) 

where we have used Tr{aia jok) ~ 2iejjK- 

Fermionic action. The coupling of Dirac fermions -0 G C°°{M,C) (E)^ to 3d gravity is given 
by the following real action 



M 



1 - - i 

- (ip eAeA \/tp - S/tp A e A e V") + 77m- Tr(e A e A e) ^tp 

2 6 



(6) 



where m G is the mass of the fermion field, the Dirac conjugate ip is given by the Hermitian 
conjugate (see the appendix) and the symbol V := V(d)) denotes the covariant derivative with 
respect to the spin connection. Equivalently, using the fact that the Clifford algebra C(3) is also 
a Lie algebra, namely ct/CTj = iej^ aK + VUj the Dirac action can be recxprcssed as 



5z3[e,w,^,V] = \j^,^^iJKe' Ae^ A Q(^a(e^^) V^- - V^a(e^)^) - ^me^ V^) ■ (7) 

The presence of the second term in the sum is to ensure the reality of the action keeping the same 
equations of motion for ip and ip. The introduction of such a term comes from the fact that in 
the presence of torsion, the usual Dirac vector current is not a total derivative [29] . 

2.2 Equations of motion 

We now consider the equations of motion coming from the variation of the total action S — 
Sgr + Sd with respect to the spin connection and to the (embedded) triad. 

The variation of the connection yields the effects of the fermionic field on the spacetime torsion 

f = -^e'jK (Tie' Ae''^^P, (8) 
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while the variation with respect to the triad encodes the effects of matter on the curvature of 
spacetime : 



From now on, we will set k = 1 and rescale the co-frame e ^ ie such that it will be regarded 
as taking value in spin(3) ~ su(2) = M.{Xj}i, with the convention Xj = iai. Having established 
the classical setting of the theory of Dirac fermions interacting with the geometry of an Einstcin- 
Cartan spacetime, we now turn to the quantisation of the fermion/gravity coupled system. 

3 Spinfoam quantisation 

Feynman. following the intuition of Dirac, realised that quantum mechanical systems could be 
approached by two different kind of formulations: the canonical or Hamiltonian approach and the 
covariant or sum over paths formalism. In the case of quantum gravity, the canonical approach 
is called loop quantum gravity and precisely predicts a discrete picture of quantum space. The 
fundamental excitations of the quantum gravitational field are one dimensional and supported 
by coloured, diffeomorphism invariant graphs. Such states are called spin networks. The world 
surface swept by a spin network forms a collection of two dimensional surfaces meeting on edges 
and vertices. The resulting object is a coloured two-complex called a spinfoam. A spinfoam can 
therefore be though of as a history of quantum space (a quantum spacetime) and a sum over 
spinfoams interpolating between two spin networks can be viewed as an implementation of the 
Misncr-Hawking three-geometry to three-geometry transition amplitude for quantum gravity. 

One of the most remarkable aspects of the spinfoam approach is the fact that a surprising num- 
ber of independent research directions converge towards the same formalism. Indeed, spinfoams 
arise for instance, independently from loop quantum gravity, as a tool of the covariant formulation 
to compute the partition function of a large class of BF-like field theories. In the case of pure 
topological BF theory, the idea is to use the lack of local degrees of freedom of the theory to 
calculate the path integral with lattice regularization techniques without needing a continuum 
limit. In the case where local degrees of freedom are present, such an approach is therefore, in 
essence, only a (non exact) discretization of the theory of interest until a sum over triangulations 
is precisely understood. A fair amount of work has been devoted to such techniques, as recalled 
in the introduction. 

Three dimensional gravity is a topological field theory and hence nicely fits into the spinfoam 
framework. In fact, the first model of quantum gravity ever written was a spinfoam quantisation of 
ricmannian three dimensional gravity: the Ponzano-Rcgge model [30j . In this paper, we generalise 
this model to include the interaction of fermions with gravity by giving a concrete meaning to 
the sum over all possible paths in configuration space weighted by the complex exponent of the 
action, i.e the path integral or vacuum to vacuum transition amplitude. 

In the presence of fermions, the theory has local degrees of freedom. Therefore the fixed 
triangulation that we will use below induces a physical cut-off of the degrees of freedom of the 
theory. In this paper, we do not address the problem of the removal of this regularization. We 
discuss some potential possibilities in the conclusion. 

This section is devoted to present the model. We first discretize the theory and regularise the 
path integral. We then integrate over the fcrmionic degrees of freedom and expand the resulting 
functional determinant in inverse fcrmion masses. This expansion yields a sum over fcrmion paths 
wrapping around the spinfoam. We finally compute the path integral of the coupled system. 




(9) 
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3.1 Discretization 

3.1.1 Manifold discretization 

First, wc pick a triangulation T of our spacctimc manifold Af [^. This means that wc subdivide M 
into a finite union of three-simphces (tetrahedra) A'^ such that any pair of such simphces are either 
disjoint or meet on a common sub-simplex (triangle or segment). Two different three-simplices 
can have at most one common triangle. All simplices and sub-simplices are oriented and the 
orientations of two triangles in the boundary of two different three-simplices along which they 
are glued together are assumed to match; the one-simplices of those two triangles are assumed 
to have same coordinate length and must have pairwise opposite orientations. If n denotes the 
number of simplices of the triangulation T, M = U^^]^A|. We will make the assumption that 
the simplicial manifold defined by T admits a spin structure. We furthermore assume that each 
tetrahedron A"^, together with the four boimdary tetrahedra glued along its four triangles, belongs 
to a given coordinate patch U C M and that all three-simplices are diffeomorphic to the standard 
tetrahedron in K'^. 

We will work with the dual complex T* of the triangulation or more precisely with its dual 
two-skeleton. This means that we place a vertex v in the center of each tetrahedron A'^ S T, 
link the vertices with edges e going through the triangles A^ € T and define a face / as a closed 
sequence of edges intersecting the segments A-'^ € T. We fix an orientation and a distinguished 
vertex for each face / and call its edges 61^,62^., ...,eAr^. taken in cyclic order starting from the 
distinguished vertex. The complex formed by the cells {v, e, f) is the dual two-skeleton of T. 

For consistency purposes that will become clear shortly, we furthermore refine the above cellular 
decomposition by introducing a subdivision of each dual face / into so-called wedges w (see |31j 
and references therein). A wedge (see figure 1) is defined by first identifying the center of a dual 
face / as the point where the latter intersects the associated segment A^ of the triangulation. By 
drawing lines connecting the center of / to the centers of the neighboring dual faces, one obtains 
a subdivision of / into wedges w consisting in the portions of the face lying between two such 
lines. In other words, a wedge belonging to a given face / and to a given three-simplex A'^ is 
the intersection between / and A'^, w = f D A'^, as pictured in figure 1. Each wedge acquires 
an orientation induced by the orientation of the face to which it belongs and is unambiguously 
associated to the simplex to which it belongs. The complex formed by these finer 2-cells, their 
boundary edges and the vertices lying at the endpoints of the boundary edges is called the derived 
complex T"*". There are two types of 1-cells in T+. We distinguish between the edges of the 
derived complex T"*" whose intersection with the edges of the dual complex T* is non empty and 
the others. We will note e+ the edges of T+ belonging to dwHdf and e_ the edges exclusively in 
dw. Note that the e+ type edges are exactly half edges of T*. There are three types of vertices 
around each wedge. The center of the face to which the wedge belongs, noted c, one vertex v 
belonging to T* and two vertices forming the intersection of the edges emerging from 1; with the 
two triangles meeting on the segment defining the center of the associated face. These two new 
types of vertices will be noted . 

We denote by V,£,T the set of vertices, edges and faces of T* while the set of vertices 

^By BF-like, we mean any theory expressable as a free topological BF part plus a polynomial function of the B 
field l§6l . 

■^Note that all compact three dimensional manifolds are triangulable. 
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V U {«+} U {c}, edges {e_} U {e+} = U f + and wedges {w} of T+ will respectively be noted 
V+, £± and W. We will call yV„ the set of wedges restricted to the wedges belonging to the 
tetrahedron dual to the vertex v. We will call n{V) = n the cardinality of the set V. 

We then assume that the global trivialising section of V is picccwise constant in each tetra- 
hedron. Hence, the spinor and vector bundles associated to V arc also globally trivialised by 
a picccwise constant section. We therefore have a basis {ei{v))i of V (or cquivalently a basis 
{Xj{v))i of spin(3)) and a basis (/q(u))q of V for each tetrahedron inside which the vector and 
spinor fields are globally defined. 

3.1.2 Simplicial field configurations 

Wc now turn towards the discretization of the field content of the theory. The main hypothesis is 
that the cellular decomposition is sufficiently refined so that the fields of the theory can be though 
of as piecewise constant inside each three-simplex A^. Let the latin indices from the middle of 
the alphabet ... and from the end of the alphabet ...,t,u,v respectively denote the vertices of 
the triangulation T and those of T+ . We denote by ab the segment linking the vertices a and h. 
For each vertex w G V, 

A simplicial soldering form configuration e is a map El 

e:Wv sptn(3) (10) 

w e,u,. 



Here, e^, := e^ denotes the line integral of e along the segment ij G T dual to the face containing 
the wedge w; eij = e(vy ) = Vij = J.^ e = Jj s*je, where Sij : J C M — > M is the segment linking i 
and j and Vij = Sij denotes the tangent vector. The components of e^, are measured with respect 
to the frame {Xi{v))i associated to the simplex containing the wedge w: Xj{v). This 

is one of the reasons why the complex T* needed to be refined to T+: a face / being shared 
by different simplices and the trivialisation of V being chosen to be piecewise constant in each 
three-simplex, the assignment of a simplicial triad configuration to the face / would have been 
meaningless: with respect to which frame would we express it ? 

The wedges of the face / are in one to one correspondence with the thrcc-simpliccs neighboring 
the tetrahedron dual to the vertex v and all share the same segment ij as it defines the center 
of the face. Hence, the vectors in spin(3) ~ M'^ assigned to each of these wedges will in fact be 
the same vector measured with respect to the frames associated to the simplices containing each 
wedge. As we will see, the quantum treatment of the theory will naturally impose on this vector 
to have same norm in different frames. 

By using the coordinates a;^ associated to the patch U Z) v , we see that, at first order in the 
segments lengths, we have ef^ ~ e^(a::«)A^- where Xy is the coordinate of the dual vertex v (it 
could be any point inside v*) and = = x^(j) — x'^(i) is the coordinate length vector of the 
segment ij. 

The spin connection uj assigns a spin rotation to the edges belonging to £"+, i.e. to the dual 
half edges (forming the intersection between the boundaries of the faces and the boundaries of the 
wedges) 

ll> : £-f- Spin{3) (11) 

e ^ Ue^. 

1/2 

The spin rotation is defined to be the parallel transport matrix Ue^ '■= Uuv{g) = uv{g) = 

^More precisely, a simplicial triad is a functor from the manifold groupoid tt(M) (the category whose objects 
are the vertices of the triangulation and the invertible morphisms are the edges linking the points) to the category 
spin(3), where there is one object and the morphisms are the additive group operations. 
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P expJ"" , g G G, of the spin connection. Wc have used the notation tt : G Aut (V) for the 
spin j representation (see the appendix), and P means path ordering. If we reverse the orientation 
of an edge, the spin matrix gets mapped to its hermitian conjugate : Ui,u = C/^^- 

To complete the picture, we also assign group elements Ue^ to the edges £- of whose 
intersection with the boundaries of the faces are empty, i.e the edges converging toward the centers 
c of the faces. These holonomies will be called auxiliary variables as they carry only information 
about the topology of the faces to which the fcrmions will turn out to be insensitive. 

By using the coordinates associated to the patch U C M containing the simplex dual to 
the vertex v, we see that, at first order in the edges lengths, the holonomy approximates the 
connection: Uuv — 1 + A(^^ct)^(a;„). 

The discrctized curvature is assigned to the wedges and defined to be the holonomy U^, of the 
spin connection along the boundaries of the wedges w : 

U^= Y[ Ue. (12) 

e^dw 

This is valid since the curvature is the obstruction to the closing of the horizontal lift of an 
infinitesimal loop. More precisely, Uw is of the form Uw = t/ei^.i7ei_C^e2-C^e2+ (see figure 1). The 
base point of the holonomy is fixed by discrctized gauge invariance arguments to be the dual 
vertex w e T* 36J. 

We now turn towards the simplicial field content of the fermionic action. Fcrmions, as sections 
of the spinor bundle, must live on 0-cells of the discrctized manifold. The structure of the classical 
Dirac action © furthermore shows us what type of 0-cells should the simplicial fields be assigned 
to. The volume, or mass, term tells us that the fcrmions live on 0-cells of T* or T"*" topologically 
dual to three dimensional regions and the area, or kinetic, term implies that the fermions travel on 
1-cells dual to objects supporting areas. By inspection, fermions can only be assigned to vertices 
of T*. Hence, 

The fermionic fields assign a spinor e V (resp. a co-spinor ip^ G V*) to each vertex v of 
V, 

^, (V^) : V ^ V (r ) (13) 

Our simplicial fermions live in the discrete analogue of the vector space of sections T{SM) ~ 
C°°(M, C) (8i V of the spinor bundle SM . This space is the vector space of V- valued sequences 
{V' I "0 • "1^ ^ V} ~ S{V) ® V mapping the space of vertices V '-^ N into the Clifford module V. 
Here, we have noted iS(V) the vector space of complex valued sequences. Wc choose as a basis 
for S{y) the set {e'"}^ of n vectors in iS(V) such that e'"{u) = (5^. Picking a basis {fa}a in V, a 
simplicial spinor field -0 G S{V) (g) V reads = i^^e" (g) /„, with 0(w) = 0^/„ G V, e C. Wc 
will note S the vector space SiV) (g) V of simplicial spinor fields. 

We furthermore need to require that our simplicial fcrmions obey the the Fermi-Dirac spin- 
statistics. For instance, the fermionic components, i.e the coefficients "0" appearing in the ex- 
pression of a simplicial spinor in our chosen basis must anticommute; tp"tp!^ + '0^''/'" =0. To 
mathematically capture this physical fact, it is necessary to change the field of numbers on which 
the vector space S is built on from the field of complex numbers C to a supercommutative ring. 
This ring is given by a Grassmann algebra that we will note Qt and which is constructed as follows 

Let (resp V„) denote the copy of the Clifford module V (rcsp of the dual complex conjugate 
Clifford module V ) associated to the vertex v of the dual complex T* . 



''In fact, the space of smooth sections of the spinor bundle r{SM) is a module over the commutative ring of 
smooth complex valued functions C°° (M, C) . 



8 



Consider the complex vector spaces E and E obtained as the direct sums 

i? = 0V,, E* = ^Y:. (14) 

vev vev 

The Grassmann algebra Qj- is obtained as the exterior algebra over the complexified realification 
of E (identifying E with its dual vector space E*): 

gr = /\(E®E*y (15) 

As a complex vector space, Qq- is of dimension 2'*" . We finally make a chart choice on Qr such 
that the algebra is defined by the generating system {V'"j V'ualf.a of the (odd) elements of Qq- 
consisting in the components of the simplicial spinors and co-spinors. 

Fermions will therefore not live in § regarded as a vector space but rather in § considered as a 
supervector space |33j . i.e a module (in fact, a bimodule) on the ring given by Qt. An element ip 
in § will be written as '0 = i^'^e" (g) fa with ^/)" G K^{E ® {0}) C Qt, where the former is included 
in the later as a vector subspace. Note that we now need to consider 5(V) entering the definition 
of § as the module of Grassmannian valued sequences. 



3.1.3 Simplicial action 

We can now construct the discretizcd version of the fermion/gravity action. As a simplicial action, 
we take the real valued functional St = Sgr,t + Sd,t, where Sgr.t and Sd,t are defined by 

SGR,r[ew, 5e] = ^Tr (e^,[/^) , (16) 

and 

\ UV V / 



where e €z £± (resp. e G £) for the gravitational part (rcsp. the fermionic part) and the sums in the 
simplicial gravitational and Dirac actions are respectively over all wedges and all edges uv linking 
the adjacent vertices u.v. In equation (jl6p . the trace is taken in fundamental representation of 
the spin group. In pT)) . 

A-uv = ^ut.(e^) = ^ ^ ei^e^'sgn{w,w') (18) 

w ,w' Duv 



with ^ 

Aiuv^-^ X! '^iJKeie^,sgn{w,w'), (19) 

w.w'Duv 

is the discretized version of the two-form e A e appearing in the kinetic term of the fermionic 
action. This term is the spin space analogue of the two-form (*e A e) where the star * denotes 
the hodge operator acting on the exterior algebra A{V) over 1/ ~ R^. At the simplicial level, 
we are looking at the integrated version of this two-form on the triangles of T or, cquivalcntly, 
we are considering the evaluation of this two form on all possible couples of vectors {vij,Vik) 
generating each infinitesimal triangle ijk of the triangulation. On a particular triangle (vij^Vik), 
{e Ae){vtj,Vtk) = {*e A e)i{vtj ,Vtk), with 

(*e A e)/(vy,-yjfc) = \eijKe^ A e^{vij,Vik) = \(-ijK{vij A VikY'^ = (wy x Vik)i, (20) 
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where (it A and X denotes the usual cross product on M'^. This evaluation therefore 

yields (the image in spin(3) of) the internal normal area one-form to the triangle ijk measured in 
the inertial frame (e/(a;„))/ associated to the simplex dual to v which is one of the two tetrahedra 
containing the triangle. In expression ([T5|) . the sum is taken over the three possible pairs (without 
counting the permutations) of wedges meeting on the edge e = uv. This explains the 1/3 factor. 
Finally, sgn{w,w') equals ±1 depending on the sign of the (coordinate) area bivector. This factor 
ensures the antisymmetry of forms at the simplicial level. 

The symbol assigned to each dual vertex denotes another polynomial function of the dis- 
cretized gravitational field which is given by 

Vy = V.„{e.uj) = Tr {e^e.u,'e.u,")sgn{w,w',w") (21) 

W.w' ,w" Dv 

= YeVs! ^ eijKeiei,e^„sgn{w,w',w"). 



This is the discretization of the integrated volume three- form Tr{e A e A e) appearing in the 
mass term of the Dirac action. Once again, sgn{w, w' ,w") = ±1 here depending on the sign of the 
volume generated by the three vectors , e^' , e^," G M'^ : it carries the sign of the orientation of 
the vectorial basis formed by the three vectors. The numerical factor comes from the fact that the 
volume of the tetrahedron dual to the vertex v appears more than once in the above expression 
Indeed, there are C| = 20 different triples of vectors constructable from the six faces dual 
to the six edges of a given tetrahedron Ej. Out of these twenty triples, only sixteen, without 
counting the permutations, do not share a dual edge, i.e span the volume of the tetrahedron dual 
to the vertex v. We call these triples admissible. For each such triple, there are 3! = 6 possible 
permutations and hence 16 x 3! times the volume of the three-simplex dual to the vertex v in 
expression (PT|) . 

The presence of this volume term in the simplicial theory has dictated the discretization 
procedure that we have followed: we have suppressed the centers c e V+ of the faces and the 
vertices of type w+ G V+ as a possible habitat for the fermions. The reason is that the centers 
are not topologically dual to a three-dimensional region and the vertices of type € V+ belong 
to a given triple of wedges which is not admissible, i.e it does not generate a three-volume. 
We accordingly see how the rigidity of the structure of the Dirac action forces us to consider a 
particular type of discretization procedure. The fermions can only he assigned to the vertices V 
of the dual complex T* and by association the spin connection can only map the edges £ of T* 
to the spin group G. 

The fermionic part (jl7p of the simplicial action is a polynomial function (in the sense of exterior 
multiplication) from the simplicial field configurations into the second exterior power 
of the Grassmann algebra Qr- We can exhibit the quadratic structure of the fermionic action via 
a compact expression by introducing a generalised Hermitian form, i.e a C/r- valued, Qq--vao(hAe 
morphism (, ) : § x § ^ Qr for which the chosen basis of S is orthonormal: 

SD,T[eu,,ge,i^v^M = {'^,W{e.^,ge)-il^) = Y^'^uWuv{e^,ge)'4'v e f\^{E®E*), (22) 

UV 

where 

W^:^p=^-^^{Ai^,(a'Uu.rp~Ai,Muo'n)-\Vu5u.5% eCcGr. (23) 
are the complex coefficients in the decomposition of the clement Sd,t hi^ the basis {ipi^a^''Pv}-u.v,ai3 

^Cn = ( ^ ) is the binomial coefficient. 



10 



of {E © E) given by the matrix elements of the t/r-module homomorphism : § ^ § on the 
space of simplicial spinors. Note that we have rescaled the fermionic fields tp — *■ y^m/3ip such that 
[iP] = m3/2. 

To prove that the simplicial action converges towards the correct continuum limit we procede 
in two steps. First, by injecting the first order expansions of the discretized fields in the discrete 
gravitational action p6p and by using the fact that the simplicial curvature associated to a given 
wedge w approximates the curvature at first order in the area of w, it is straight forward to see 
that Sgr,t —Sgb. pointwise when n — > oo and the size of each tetrahedron shrinks to zero 
uniformly in the formal symbol k labelling the three-simplices. 

Secondly, we analyse the fermionic sector Sd,t- Once again, we expand the simplicial fields in 
small vertice to vertice inter-distances corresponding to a high degree of refinement (n >> 1) of 
the triangulation T. The Taylor expansion of the fermionic fields around the vertex u £ V yields, 
in the coordinates associated to the patch U C M containing the vertex u together with its four 
boundary tetrahedra, 

V'. =V'« + A(J„9^Vu+0(Aa;2). (24) 

By using Stokes theorem, it is furthermore possible to show that for a fixed vertex u e V, the sum 
of the four coordinate area bivectors corresponding to the four triangles dual to the four edges 
emerging from u sum to zero |37| . This is due to the fact that the integral of the unit normal 
pointing outwards to a two-sphere piecewise linearly embedded in R'^ must vanish. Explicitly, 

VueV, E («^^- ^ "»'^)^'I-=«.* = 0' (25) 

V 

where ijk is the positively oriented triangle dual to the half-edge uv Cz 

Finally, we use the fact that the normal coordinate area one-form to a given triangle of T 
and the dual length vector going through generate a coordinate three-volume [31] : 

Y,A^uvAx'^,=2n{M)S;, (26) 

uv 

where n{M) is the coordinate volume of the spacetime manifold M and the normal one form 
Afiuv is related to the area bivector {vij A VikY" of the triangle ijk dual to the edge uv through 
a (coordinate) hodge operation: [vij A VikY" = ef^^PApuv 

By inserting the Taylor expanded fields in the discretized action (fT7|) . using expression ((25|) 
to kill the undesired terms and by finally recovering the coordinate volume trough ([26]) . we can 
conclude that, if the sequence of simplicial manifolds considered here converges to a smooth 
manifold with spin structure, the simplicial action (jl7p converges pointwise towards (minus) the 
Einstein/Dirac continuous action when the characteristic length of the simplices gos to zero and 
the number of simplices tends to infinity, at least formally. 

We finally verify that the action (fT7|) is real as it is the case in the continuum. 

3.2 Fermionic integration 

We are now in position to give sense to the formula ([1]) by using the discretization prescription 
discussed above. The main obstacle to analytically work with path integrals resides in the defini- 
tion of the functional integration measure. The lattice regularization used here provides us with 
an effective tool towards such a definition. 
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3.2.1 Berezin integral and fermionic determinant 

The formal expression ([1]) becomes: 




dfxi^,, V.)) e*^-[-""3-'^-'""l, (27) 



where de^ is the Lebesgue measure on g = spin(3) ~ M.^, dge {due) is the normalised left Haar 
measure on G = Spin(3) (which coincides with the right Haar measure because G is compact: in 
this sense, G is said to be unimodular) and the symbol d^[tp^, ipi,) denotes the Berezin integral 
[55] on Qr- The Berezin integral is the element of the algebraic dual of the Grassmann algebra 
Gt defined by 



Vae^r, / «rfM(^., ^.)-(riri^) (nn^lr) « eC. (28) 

"'St \„=„ a=2 '-'Vva / \v=n a=1 / 

As a direct consequence of the above formula, the Berezin integral of the exponential of a quadratic 
expression in the Grassmann variables yields the (totally antisymmetric) coefficient in front of the 
basis element J^^ ^ ip" Yiv a ''Pva the (finite) development of the exponential. This means that 
the Berezin integral in ([^7]) produces a sum over all possible alternating 2n-tli powers of the 
coefficient W, i.e the determinant of W regarded as an endomorphism of the t/r-module §. This 
exact manipulation is due to the fact that we have cut down the number of degrees of freedom of 
the theory from infinite to finite by introducing a lattice discretization. The computation of the 
Berezin integral appearing in the partition function follows 

/ dM^,, ^„)e'(^'^^) =detW^, (29) 

and we can accordingly integrate out the fermionic degrees of freedom of the theory hence obtaining 
the usual path integral for gravity with insertion of an observable corresponding to the functional 
determinant det W: 



Z{T) = <detW >T (30) 
= (llj^de^^ {llj^du,_^ (llj^'^de^^ [detW{e^„ge)]e'^-^'-^--''-\ 

We now develop a tool to calculate the functional determinant det W leading to the computation 
of the path integral of the model. 



3.2.2 Fermionic loop expansion 

In statistical physics, the high-tcmpcrature expansion is a valuable tool for investigating critical 
phenomena. In lattice-like theories incorporating massive fermions, an analogous procedure called 
the hopping parameter expansion [39j exists. The technique consists in expanding in inverse 
fermionic masses in the heavy mass regime of the theory. Of course, the sole consideration of the 
lowest orders in the expansion is reliable only if the fermions are heavier than the energy scale 
of the other fields of the theory under study. In the case of quantum gravity, the typical energy 
scale is given by the Planck mass which implies that the expansion parameter turns out to be 
large. However, if the expansion is finite and considered as a computational tool and not as a 
perturbation theory, we believe that the hopping parameter expansion in quantum gravity remains 
a valuable tool. This idea has in fact successfully been implemented in |40j where the coupling 
of massive fermions to two dimensional Lorentzian dynamical triangulations is considered. The 
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large value of the expansion parameter simply implies that only a complete consideration of all 
orders will yield a meaningfuU result. In other words, the hopping parameter expansion can 
only be used in quantum gravity as a computational technique, not as an effective description 
of the theory: we will not construct a perturbation theory but use the finite hopping parameter 
expansion to compute explicitly the partition function of the theory. The crystal clear geometrical 
interpretation of the hopping parameter expansion in terms of fcrmionic loops where the fcrmionic 
couplings to gravity appear as extremely natural and the explicit solvability of the model comfort 
us in the idea that this technique is an appealing tool to couple matter to gravity in the quantum 
regime, at least on a fixed triangulation. Indeed, it is important to stress that the finiteness 
of the expansion that we are about to derive is due to the regularization defined by the fixed 
triangulation. The behaviour of the expansion under refinement of the triangulation is a question 
that we will leave open. Leaving this issue aside, we now show how to compute the partition 
function of the theory order by order in the expansion parameter obtaining a reliable evaluation 
of the fcrmionic determinant. 

The functional determinant that we wish to calculate is defined by equation (j29p . The first 
step is to realise a crucial symmetry of the W endomorphism under charge conjugation. The 
charge conjugate spinors 

^^ = C-V* and V^^ = ^*C, (31) 

where C is the quaternionic structure defined in the appendix and * denotes complex conjugation, 
are solutions to the Dirac equation with charge sign inversion e — > — e (when coupled to an 
electromagnetic field). 

The Clifford representation p : C{3) ^ End{€?) is equivalent to the Hermitian conjugate 
representation acting on C^* via the bijective intertwining operator!^ 

A = €oC:<C^^C^\ (32) 

where e : — > C^* is the two dimensional symplectic form (e : C'^ ^ C^* is the complex conjugate 
map) whose matrix representation is given by the totally antisymmetric tensor 

This vector space isomorphism mapping onto its contragredient representation space realizes 
the equivalence between the two representations. Since A is the identity map, ^(7)^ = a"^ = a, we 
find that the charge conjugation matrix C is given by the inverse of the dual pairing map e. 
Putting indices, the conjugate spinors and co-spinors are defined by 

Ve"=e"''V^^ and = -^-^a, (34) 

where ei2 = e^^ = 1 and Cape^'^ = —S'^. One can verify that three dimensional Euclidean charge 
conjugation is indeed an anti-involution, preventing us from defining true Majorana spinors. 

These results trivially generalise to the case at hand namely the simplicial, spacetime depen- 
dent, Grassmann valued case. However, C now becomes a Qx-module morphism C : S ^ S, where 
the conjugate module § is defined with an appropriate real structure on the complex vector space 
Gt- 

Remarkably, the fcrmionic endomorphism W transforms as 

CWlC-'=W,,.u ^eWuve^-W^^, (35) 

"An intertwining operator if> between two representations Try and ttt^ of G acting respectively on the vector 
spaces V and W is a, linear map between V and W satisfying <f> o Try = ttw ° The space of intertwiners is a 
vector space noted HomciV, W). 
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under charge conjugation. 

In order to exploit this symmetry, we perform an affine transformation on and introduce a 
family (jji, i = 1,2, of symplcctic Majorana spinors (rcsp co-spinors) : 

0ic = <^ij4>j, (36) 



defined by 



(37) 



62 = -K^/'c-iV'), ^2 = ^(V'c + ^V') 



The components of these new fields are not independent. The symplectic Majorana condition 
implies the following set of constraints on the field components 



"2 - Vip 



IZ (38) 



Hence, we have not modified the number of degrees of freedom involved in the theory. As we are 
about to see, the introduction of the symplcctic Majorana fields however simplifies the situation 
rather drastically. 

Indeed, reexpressing the Dirac action in terms of the new fields, using the conjugation sym- 
metry psp . the above constraints and the Grassmann nature of the field components, we see that 
the two families decouple 

SD,r['>IJv,i'v\ = '^i'^Wuvi^v = -«(5r[0i«] + '5'r[02i.]), (39) 

where 

St[(I)v\ ■■= ST[ef,ge,(j)v] =Y,KWuvc.i3iew,ge)(l)^, (40) 

uv 

and Wuvafi = ^a-yW^^ p. This action is very similar to the original simplicial Dirac action except 
that the degrees of freedom are reduced by a factor two via the symplcctic Majorana constraints. 
The sum of the two actions however restores the original degrees of freedom, they have simply 
decoupled. This procedure gives us the opportunity of calculating the fermionic determinant 
exactly. With the original framework in terms of Dirac spinors, the presence of "too many" 
degrees of freedom per spacetime point prevents us from computing the terms at all order in the 
expansion that we now define. 

By writing the integral definition of the fermionic determinant (|29p in terms of the new fields, 
we obtain 

det W = Pf{eWf, (41) 



where 

Pf{^W) = ^ f dA^(0„)e^-['=™^^-^"], (42) 
2" Ja{e) 

is the Pfaffian[3 of the antisymmetric endomorphism eW, Wuvaf3 — ^WyufJa, whose square yields 
the determinant of eW equating det W because the determinant is a morphism of endomorphism 
algebras. Note that the Grassmann algebra integrated over is now "half" the original Grassmann 

^The Pfaffian of an antisymmetric matrix M^u of rank 2n is defined as follows. Let M = 5Za<i) ^'^ab £a A e;,, 
where {ea}a—l,...,2n denotes a basis of M^", be the associated bivector. The Pfaflian Pf{M) is then defined by the 
equation : i^M'^'^ = Pf{M) ei A 62 A ■ ■ ■ A 62,1. 
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algebra C/t, namely the exterior algebra over E. Furthermore, we have chosen as a basis for E 
the components of the symplectic Majorana field (j>. 

We now first concentrate on the calculation of the Pfaffian of the endomorphism eW before 
computing the square. The idea is to construct a finite perturbative definition of the integral 
in terms of inverse fermionic masses. We rewrite the Pfaffian by splitting the exponent into the 
volume contribution and the kinetic term of the W endomorphism (see ()23p ). The last is then 
expanded in inverse fermionic masses. We obtain the finite expansion 

^/(e^) = ^E«'rfc, (43) 

where the coefficients Tk := Tk{e.w, ge) are given by 

^fc 77 XI / dfi{(l)i,)T\e~^^'''^''^'^'' [(j)u,£Du^y^(l)y^...(j)uf,eDu^y^(j)yJ , (44) 

i^- ..... .. „ A(ij) V 



tll-Ul ,...,tlfct)fc 



with the kinetic term D 



3i 

Duv = {AiuvCt'Uuv ~ AiyuUl^a') . (45) 



Putting back physical units (c is left fixed to unity), the parameter a is given by a = = ^ 

with Ip = Gh and rrip = ^ respectively denoting the Planck length and mass. To explicitly realise 
the appearance of the Planck length, we have rescalled the triad such that e Finally, the 

passage from the exponential of a sum to the product of exponentials is due to the sole presence 
of even elements of A {E) in the exponent. 

We are now ready to perform the integrals order by order in a. Most terms in the expansion 
vanish after integration. The Bcrezin integral yields a non- vanishing result if and only if each vertex 
w £ V is exactly saturated by two (non identical) symplectic Majorana fermionic components. As 
a result, each vertex either carries a volume contribution or must be a source and a target point 
for a fermionic worldline. For instance, we obtain the following identity 

/ rfM0.)n^"*''"*"'*"^*>^=(-i)""'n^-'^*«^"''- (46) 

Jme) 



The generalisation of the above formula to an arbitrary number 2k of field insertions implies that, 
leaving aside the sum over vertices, the resulting term in the Bcrezin integral of the coefficient 
of order k ()44p is a contribution consisting in a closed sequence of products of kinetic terms D 
going through k vertices of the dual complex T* together with volume weights living everywhere 
in the complementary of the sequence. These loops can close on a vertex u and reopen on the 
neighbooring vertex v yielding disconnected sequences such that the number of vertices traversed 
by the connected components sum up to /c. A sum over the k vertices supporting the closed 
product is then implemented by the sum over vertices in equation (j44p . 

Forall k e]0, n[ and for any fixed set Vk = {f i, W2, ffe} C V oi k two by two adjacent, distinct 
vertices, such a contribution, noted y^j^^'"''^''\ will be defined in the connected case by 

ri"^-""''-' - (-l)'=-^(-l)"-''^ n TrD{D.,..D..,,,...D,,.,). (47) 

where Tr^j : End{S(y)®Y) End{S {V)) is the partial trace acting on the Dirac part of End{S). 
If the sequence Y^f^^''"'^''^ is not connected, there will be one Dirac trace appearing in the above 
formula for each connected component. 

The coefficient F^ is accordingly a sum over connected and/or disconnected contributions : 
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(a) 



(b) 



Figure 2: (a) A fermionic loop of order six. (b) An order twelve product of two loops. 



r^-^[EEri"'-"^\ (48) 

r 

where the first sum is over possible vertices supporting the loop and the second formal sum symbol 
implies a sum over all possible connected and/or disconnected contributions saturating the k 
vertices Vfc. 

We now compute the extremal configurations. The order zero is the no loop case and is given 

by 

ro = (-i)"n^- (49) 

V 

Note that the first order yields the trace of the D endomorphism by virtue of equation (1^5]) which 
is null. 

If the order of the loop k equates the number of vertices n, there is no place for the volume 
contributions and the order k coefficient in the expansion (j43p yields 

r„ = T{-iy'-'TrD (-D„,.,-D„,.3...i?.„.J . (50) 

Once again, the loop can be disconnected and built out of lower order connected components in 
which case there is one Dirac trace per connected piece. What do these loops have to do with 
vacuum fermionic ? Each Dirac trace over a closed sequence of D endomorphisms generates 

a sum over 2^ — X)p=o different contributions by virtue of equation P5|) . For a fixed set 
vi,V2, ■■■,Vk of k vertices, 

TrniD,,.,D,,.,...D,,.,) = |^ (51) 

i—l 

Each element 7 in the above sum is a closed fermionic path going trough the k vertices vi, 
Such a loop is attached to a decoration, i.e a pattern of sigma matrices insertions together with a 
polynomial function of the discretized triad Cw 

Let us introduce some terminology. As suggested by the notation above, we note 7^ a loop 
wrapping around k vertices of T*. Such a loop will be called of order k. All fermionic loops 

®More precisely, by denoting 7v^, the Dirac trace Tr£i(_D„j^„2 _D„2U3 •■■^Ufcui )i could define the cocfiicient Fj. 
by introducing a partition lip of Vfe with p subsets: Vfe = U^_-^Vk., — which case 

where the first sum is over all vertices, the second over all numbers of subsets p constituting the partition and 
the last sum is over all possible partitions with p subsets. 
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are connected, do not go trough the same vertex more than once and when a product of loops 
is generated at given order in the expansion by a contribution pj^"!' - '"'--)^ ^j-^q different loops do 
not touch each other (see fig.l). For each fixed set of k vertices vi, ...jffc € V, the contributions 
appearing in the above sum (|5ip differ between themselves by their decorations, i.e by the posi- 
tioning of the sigma matrices inside the trace and their pattern of contraction with the (internal) 
area co- vectors A/. We call such a contribution 7' in the above sum a decorated loop and note 
£fe the set of decorated /e-loops. We define an equivalence relation ~ on the set >Cfe of loop config- 
urations of order k by calling two contributions equivalent if they can be mapped onto each other 
by relabelling of the vertices. This is the case when the two decorated loops traversing the same 
set of vertices are built out of the same sequence of patterns alternating the sigma matrices and 
the group elements inside the Dirac trace up to cyclicity. 

A typical generic decorated loop of order k looks like this: Vvi 7^ Vj, 

7fc ^ ~ ^-fi viV2-^i2 V2V3 ■■■-''^ik Ufc-oi Tru (^(7 ^Uy-^y^^ ~Uv2V3---'^ ''Uvi.vi') ■ (^2) 

This term is the first in the sum over decorations ([5T|) and will be the generic loop with which we 
will work to build the model. All of the other decorated loops of a given order k can be obtained 
and usefully classified trough the following prescription. 

The traces!^ are built out of a succession of two types of patterns of the form aU and C/V 
(see We call the first pattern (resp second pattern) a pattern of type a (resp of type b). A 

contribution of order k will contain a succession of A: — m (resp m) type a (resp type b) patterns. 
There are C™ possible combinations of such type which can be classified as follows. First, represent 
each pattern by a labelled point on a circle (type a or type b) and each combination of such patterns 
(up to cyclic permutation) by a set of labelled points on a circle. Changing the starting point on 
the circle corresponds to a relabelling of the vertices. Hence, each circle corresponds to an element 
in Ck/ ~. Our classification is established by identifying the finite number TV of possible combined 
patterns of the form "type b times type a" , i.e U^acrU, constructable out of the fixed number k~m 
(resp m) of type a (resp type b) patterns. For each fixed there are N sequences of the form 
aUaU...aUW aU^ a...W a constituting each circle. The different possible sequences for fixed N 
correspond to the different ^^-inequivalent contributions containing N patterns of two adjacent 
Pauli matrices inside the Dirac trace. We will call the sub-sequences of the form aUaU...<jU (resp 
C/V[/V...C/^cr) type a (resp type b) sub-sequences. 

The contributions to the Pfaffian have now been defined and classified. To summarise, 
the calculation of the Pfaffian yields a sum over products of connected decorated loops wrapping 
around more and more vertices of the dual triangulation T* with the increasing of the order of 
the expansion. This sum stops with an extremal decorated loop saturating all the vertices of T* . 
At given order fc, the coefficients in the expansion are symbolically written 

Ffe = ^c^, Jl K 7fc, c^fc e C, (53) 

where the sum is over all possible collections 7fe := U,7i of decorated loops 7 contributing to the 
order k. Our conventions are such that 70 = 70 = 1. We call the ensemble 7^ a loop configuration. 
The complex number c^^ contains all the numerical factors including the multiplicity of the loop. 

We can now compute the square of the Pfaffian to obtain the functional determinant det W 
and restore all the degrees of freedom of the (discrctized) theory. The determinant expansion will 
accordingly involve a sum over quadratic terms in the loop configurations. At a given order fc, 
leaving aside the sum over vertices, we obtain terms of the form ^p+q = F^.F^ with p + q ~ k. 
We distinguish between three different types of contributions. 

First, there are the terms where both p and q are smaller that the number n of vertices in 
which case the product of loop configurations generates configurations where all vertices in the 
complementary of the union of the two fermionic loop configurations are associated to square 
volumes while the vertices traversed by the fermionic loops carry a single volume contribution. 

®We are leaving the area co- vectors out of the discussion for the sake of simpUcity. 



17 



We can now find vertices where more than one fermionic Une converges and edges supporting two 
fermionic hnes. These particular vertices will be called degenerate and do not carry a volume 
weight. If p is smaller that n and q equates n, the intersection between the two fermionic loop 
configurations is necessarily non empty and the resulting loop configuration is degenerate. All 
non degenerate vertices contribute a volume term while the degenerate ones will carry a trivial 
weight. Finally, if both p and q are equal to n, we obtain a contribution consisting in no volume 
term and a double fermionic loop going around all of the two complex. 

The structure of the expansion being established, we can rewrite the partition function of the 
simplicial theory as the following finite sum 

2n 

2(^) = ^E"'^(rfe,r), (54) 

fe=0 

Here, Tk — J2p q ^p^q^ such that p + q = k, and Vfc e [0, 2n], 



(55) 



7fe 



7fc 



'7fc 



Mh) (n//--) (pX"" 



Y[fv,jkiew)'^k{ew,ge) 
L V 



where 7fc = 7p U 7^ with p + q ^ k and fv,^^ is a polynomial function of the triad whose exact 
nature (identity, volume, squared volume) depends on the vertex v and on the type of decorated 
loop product 7fe according to the discussion above. 



3.3 Fermions and quantum gravity 

We are now ready to perform the integrations in the path integral and solve the model at the 
quantum mechanical level. It is rather straight forward to compute the integration over the 
discretized connection, i.c over the group elements attached to the edges. The difficulty resides in 
the calculation of the integral over the simplicial triad field. The idea that wc will follow, due to 
Freidel and Krasnov [35], is to use the fact that any polynomial function of the discretized triad 
appearing in the partition function can be obtained by source derivation on a given generating 
functional. 



3.3.1 Generating functional 

We firstly study a modelling case where the loop configuration 7^ is of the form 7^ U 70 where 
7i; := ^("i' - '^'')^ Vfc < is the single, connected decorated fc-loop defined by (f5^ . This prototype 
situation retains all the important features of the model and will lead us to derive the full model 
containing all the other types of loop configurations at the end of this paper. Troughout this 
section, we will assume that k > 2. 

Regarding the framework developed above, the associated weight in the partition function is 
given by 



(56) 



n (^-) n ^-(^-) n ^^e,e(e.„) ( n '^'°ue) 



"67* 



ee7fc 
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The technology of generating functionals gives us the possibiUty to compute such an expression 
needed for the construction of the model. 

The idea is to rewrite the above quantity of interest as 



v^lk '"£71- ee7fc 



(57) 



j=o 



Here, a degree n element O of the algebra of monomial functions over spin(3)^" is transformed 
into a source differential operator of order n O := 0{jj) as follows 



0(e„) 



O 



i 6 



(58) 



2 5Ji, 



The generating functional Z^^\j,^k,T) is a map from the space of spin(3)-valued two forms into 
V^^ defined by 



z(^)(j,7fe,r) = (n^'^^-) (n/^^"-j (n^^ff-^ 



piS'r,j[e™,ge,e-'" 



with (/) := (/ei,/e2: ■•■1^6^) denoting a cumulative index and 



(59) 



(60) 



where '■= J, with J = J^Xj, being a g-valued current two-form. Using the traceless property 
of the elements of spin(3), we see that St,j is a discretized version of the continuum source action 



and that the operation 



6Ji 



Sj[e,uj,J]= Tr {eAR[u])+Tr (eAJ), 

J M 



(61) 



(62) 



approximates the continuum operation of source derivation at first order in the segments G T 
lengths {Uw ^ 1). This motivates our definition of the expectation value of any degree n monomial 
in the simplicial triad : 



< O(e^) >r (-)" 



S 



SJi 



(63) 



j=o 



We now calculate the above generating functional ((S^ to be able to make sense of equation 
following the next lines . 
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Simplicial triad integration. First, the Lebesgue integrations over M'^ ~ g can be performed 
yielding one delta function on the spin group for each wedge w. Then, using the fact that the 

^ j 

characters Xj '■ G ^ C; g ^ Tr n (g), provide a vectorial basis for the space of central functions 
on G and that the delta function is central as a distribution, wc obtain 



n//"-) (n 



e67fc 



J W j„eN/2 



(64) 



where dim{jw) — 2ju, + 1 is the dimension of the representation jy^ assigned to the wedge w. 



Integration over the discretized connections. The next step consits in performing the Haar 
integrals. All 1-cells belong to more than one 2-ccll. The e_ type edges are shared by two wedges, 
the e+ type belong to three or, as we shall see, four wedges or faces. The group elements assigned 
to the edges will therefore always appear in more than one trace in We accordingly need to 
develop the necessary technique for integrating over G the tensor product of two, three or four 
representations of a given g € G. 



Forall triple of unitary, irreducible representations tt, tt, n of G, let 

k i j 

:V^V (8) V denote the Clebsh-Gordan intertwining operators 



:V ® V- 



and 



These maps are uniquely 
defined up to normalisation since the vector space of thrce-valent G-intertwiners is of dimension 
one. The phase and the sign of the Clebsch-Gordan maps is fixed by requiring their reality, asking 
= and by following Wigner's convention [42]. Wc can define a Hermitian form <, > 
on the vector space of intcrtwiners with respect to which the Clcbsh-Gordans are normalised to 
unity K"^,^ >= Tr(1'«'t) = Tr(^'$) = 1. 

We introduce the non canonical bijective intertwining operator 



: V 

j 



3 

V* 

,i \ i b 
ej(ea) = e ejb 



(65) 



with ejab = (—I)-' '^Sb.^a and ei/2 := £• This vector space isomorphism defines a scalar product 
on V via (u, w) = (e(w))(w), Vw,w GV, and proves that any representation j is equivalent to its 
contragredient: ^ * = n = ej tt ej^ . Using this isomorphism, it is possible to raise and lower 

i j k 

intertwiner indices i.e to construct, for instance, elements of Hom^iy ® V, V*) : 



^^k^^ = Efc o * 



(66) 



^''More precisely, these integrals produce delta functions on G = SO(3) 1181 . In fact, we are here considering 
integrals of the form 

H 

where e{g) is the sign of cos 6 in the parametrization g{9, n) = cos 8 1 + i sind n.a of the spin group, as a mean to 
construct delta functions on G = Spin (3). 
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We will extensively use a symmetrized version of the Clebsch-Gordan intertwiner called the 
Wigner three-j map 

Lijk : V ® V «) C. (67) 

Its relationship to the Clebsch-Gordan maps is given by the following evaluation called a three-j 

I I 

symbol. Using the basis {ea}a of V, with I ~ k, 

irjki^a ®ih^ec)^(^l I ^ ^ = e"^*-^-'^) (e, , ^%iea ® St)) G M- (68) 

This intertwiner satisfies, among other properties, the symmetry conditions Lkij = ^jk ~ I'jki — 
^nT{i+j+k) ^^^^ and is normalised: < l,l >~ 1. 

We are now ready to integrate out the group elements. The first step consists in taking care 
of the auxiliary variables u. This is achieved by remarking that this operation will require the 
integration over the tensor product of two representation matrices since the internal edges to 
which the auxiliary variables are attached are shared by precisely two wedges. 

i j ■ ■ k 

Using the complete reducibility n (g) tt = tt ^'^-j of the tensor product of two irreducible 
representations, it is immediate to show that 



dg^(5)®^(ff) = n*%-: (69) 

G 



where vE'°j(e (g) e) = — ^====(ei(e))(e) — — ^==— (e,e) if i and j* are equivalent and yields zero 

otherwise. The fact that the orientation of the wedges is induced by the orientation of the faces 
to which they belong implies that the two wedges meeting on a given edge e_ will always have 
opposite orientations. Accordingly, the integration over the auxiliary variables will impose each 
wedge belonging to a given face / to be coloured by the same representation jf. Recalling that the 
representation associated to a wedge is interpreted as the quantum length number of the segment 
defining the center of the face containing the wedge as measured from the frame associated to the 
wedge, we obtain that the vector associated to the segment has the same quantum length in all 
frames [5^ . 

The generating functional ([59]) consequently yields 



(70) 



c G 



Tro ill cr'^Ue) x 



where the group clement insertions are along all the vertices v £ V oi a, given face / (m denotes 
the number on vertices, i.e the number of wedges contained in /) and xif) denotes the topological 
invariant Euler characteristic of the face. 

The second step consists in Haar integrating over the group variables assigned to each edge 
e of T*. The dimensionality of our spacetime implies that each edge e is shared by exactly 
three geometrical faces (each triangle is built out of three segments). However, if a fermion is 
travelling along the edge e, there will be an extra face, that we will call virtual, in the fundamental 
j = Jo = 1/2 representation corresponding to (a piece of) the fermionic loop. We now distinguish 
between the two cases : the purely gravitational situation, away from the fermionic loop 7^ , and 
the edges supporting a fermion path. 

If no fermions are present, each group element will appear three times in the product over the 
faces in (j70p . If i, j and k denote the three geometrical colourings of the three positively oriented 
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faces meeting on e, i.e the three lengths quantum numbers of the three segments building up the 
triangle dual to the edge e, the integral of the tensor product creates a three-j symbols together 
with its adjoint intcrtwincr 




dg TT (g)(g) TT (g)(g) TT {g) 




(71) 



where l := t^fc. If one or two out of the three faces are not positively oriented, one uses the 
isomorphism (|65p on both sides of the above equation in order to obtain the integral of interest. 
The raising and lowering of indices on the right hand side yields the correct intertwiner via (|66p. 

Let us now turn to the case where a fermion is travelling through the triangle dual to the edge 
e, i.e e G 7fc. There are now four representations associated to three geometric colourings of the 
three faces meeting on e plus a matter colouring of the virtual face defined by the fermionic loop. 

We are here contemplating the creation of quantum torsion by the fermionic field dislocating 
the quantum mechanical spacetime in which it is evolving. Indeed, if we consider the triangle 

S T traversed by the fermion line as a physical system determined by the (physical) lengths 
of its boundary segments, we can apply to it a Kirillov-Kostant geometric quantisation procedure 
obtaining the so-called quantum triangle [43| . The lengths of the boundary segments of the quan- 
tum triangle are quantised and are encoded in the associated representations. When no spacetime 
torsion is present, the quantum version of the Gauss law tells us that the three quantum length 
vectors must close. In presence of fermions, the Gauss law acquires a source term. Heuristically, 
we see that the presence of the fourth virtual face is a gap in the closing of the quantum triangle 
whose magnitude is given by the associated representation, i.e jo = 1/2 in Planck units. The 
inspection of the classical equations of motion shows that this gap is the quantum analogue of the 
failure to the coplanarity of the image under the gravitational field e in incrtial space of the three 
(coordinate) length vectors of the triangle traversed by the fermion. 

In addition, equation (|59p teaches us that, for the generic case that we are considering here, 
each fermion path along an edge comes accompanied by a sigma matrix. Recall that, classically, 
the sigma matrices are matrix representations of the embedding of the co-inertial space V* in the 
Chfford algebra C(3) and hence solder between the spin space V, i.e the j = 1/2 representation 
space of G, and the inertial space V, i.e the j = 1 representation, created by the gravitational 
field at each point of spacetime (see the appendix). In the quantum framework developed here, 
we see that these representations of C(3) are attached to the edges of the fermionic loops. How 
do they relate the spin space to spacetime in the quantum regime ? The answer to this question 
resides in the rethinking of the sigma matrices as elements of Hom^{V* ,Y (8) V*) by virtue of 

the defining equation (|109p of the spin group G. Because V* , V and V are isomorphic as vector 
spaces, we can conclude that the sigma matrices are intertwining operators between the adjoint 
representation and the tensor product of a fundamental and its dual representation. 

More precisely, by using the isomorphism Hom^iV* V*,V*) ~ Houiq^V* , End(V*)) given 
by $(e^ (8) e") = -^(o'(e^))*(e"), we obtain the following relation between sigma matrices and 
three-j symbols 



where we have omitted the representation labels of the intertwiners as there is no ambiguity and 
the scalar product used is the one induced on V* par the previously defined scalar product on V. 
This reformulation will represent the quantum analogue of equation ()109p . 

Hence, when a fermion line jo = 1/2 comes attached to a sigma matrix, we will have to 
repeatidly use the following relation 




a 



1/2 




(72) 




{g)(E) n {g)(E) n (g)® ^ (g) = ^ (a^ o 4) l. 



S 



(73) 



where :— a{e^) 1 1 1 and the sum on the right hand side is over the representations 



22 



Figure 3: Graphical picture of the emergence of the fermionic seven-j symbol. 



appearing in a chosen decomposition of the space of quadrilinear invariants. For instance, if we 

choose Hom^iY V V «) V, C) = 0^ Hom^{Y Hom^iY V, V*), the four-valent 

intertwiner appearing in the above equation is given by : 



sjk- 



(74) 



The coefficient appearing in the change of basis between the chosen decomposition and the 
one given by coupling the representations jo ^ 3 and z — fc is given by a six-j symbol. 

Integrating out all group elements and putting everything together, we obtain the following 
generating functional 



xif) 



^ if 



n{6jK.(j)- 



(75) 



Here, {6j}u(J) is the source six-j symbol defined by calculating the usual {6j} symbol out of 
the six representations associated to the six wedges sharing the vertex v with the insertion of the 
group elements e"^ . For a given vertex v belonging to six wedges coloured by the spins ji, jei 
with associated sources Ji, Jg and for a given orientation configuration of the wedges ; 



{6j}(J) = { f f f (J) (76) 

Ji J5 36 ) 

ai 02 «3 W jl 05 je \ f ji J2 "^6 W 04 js .?3 

jl 32 J3 / V a'l is a'f^ ) \ a'^ a'^ je J \ ji ^ a'^ 



The six-j symbol (with source insertions) can be conveniently rewritten as a graph reproducing 
the tensor contraction pattern. This graph is built out of four nodes n and six links I respectively 
coloured by intertwiners, issued from the Haar integrations, and irreducible representations of G 
coming from the colouring of the faces meeting on the vertex v. The resulting tetrahedral graph 
is dual to the original tetrahedron of the triangulation T; his nodes (resp links and triangles) 
correspond to the triangles (resp segments and vertices) of the original tetrahedron. This means 
that the three links emerging from a node of the spin network correspond to the three segments 
building up the corresponding triangle of the original tetrahedron. Such a graph defined up to 
topology and encoding combinatorial representation theory data is called a spin network and can 
be pictured as follows 
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{6j}(J)= ^ 



(77) 



We will fix are conventions such that all spin networks will be built out of normalised intertwiners. 

By {7j{'^)}l{J), we mean the modified {6j} symbol calculated by taking into account the 
fermionic line (and its associated sigma matrix) going trough the vertex of interest. Figure 3 
pictures the integration over the group variables yielding the fermionic seven-j. For a given vertex 
V traversed by a fermionic line and belonging to six wedges coloured by the spins ji, je with 
associated sources Ji, ...,Je, the fermionic seven-j yields 



fli 02 03 

ji 32 h 

ja ji .h ae 

a 04 a'2 je 



Ji 05 JG 



Or 



J5 

(3 04 js 
jo j4 05 



(78) 



J3 



1 /Jo i j k 

where tjgijk [ ea® e.a® e.h® Ec) 



Jo 
a 



Equivalently, the scvcn-j is defined graphically by : 




(79) 



If fc = 2, the one half line in the fermionic scvcn-j picture closes by going from one node to itself. 

It is interesting to stress that evaluating the obtained partition function at J = yields the 
expectation value of a fermionic Wilson line (with sigma matrices insertions) coupled to three 
dimensional quantum gravity. 

3.3.2 Grasping operators 

Recall that for a given fermionic loop union 7^ — 7fc U 70 of order fc, with jk denoting our generic 
decorated loop defined above, the weight in the total partition function ([54|) is provided by ([57|) . 
with a by now precise definition of Z'^-* ( J, 7fc, T). The last step hence consists in understanding 
the action of the differential operators acting on the generating functional. The key point is 
the following. If J is a g-valued current, J = Xj, then the action of source derivation on an 
exponentiated current in the representation j yields 



(80) 



i.e the image of the I-th generator of a given (dual) basis of g under the representation j of the Lie 

j . ^ 

algebra. The crucial point relies on the fact that the Lie algebra homomorphism tt, : g ^ V (g) V* 
is an intertwining operator between the adjoint representation, a j representation and its dual. 
We note 3> the associated normalised intertwiner; 



(81) 



24 



where OUf := 6(1, j, j)^ 4,, ^,>^ j{j + l){2j + 1). 

Therefore the action of source derivation with respect to the I-th component of the current, also 
called grasping, modifies a spin network with source insertion by creating a new vertex attached to 
an open line coloured by the index / in the adjoint representation. From now on, we will represent 
a link in the adjoint representation by a dashed line: 




The no fermion case. For all fermionic loop ensembles jk of the form 7p U 7^, where p and q 
are smaller than n, there is a squared volume living on each vertex where no fermions are present. 
This is a remnant of the presence of the mass term in the Dirac action. The squared volume of 
the three-simplex associated to a vertex v £ T* is calculated by acting twice with the volume 
operator 



8.16.3! 



E 



w.iu' ,w"Dv 



5J, 



w"K 



-sgn(w, w' , w"), 



(83) 



on the six-j symbol with source insertions {6j}u( J) ([77|) living on v and by evaluating the result at 
J = 0. The action of Vy associated to a given triple of wedges will create one intertwiner on each 
of the links I, I' , I" corresponding to the wedges w, w' , w" mapping the dual adjoint representation 
space into the second tensor power of the representation space associated to the link. The three 
adjoint representations spaces then get mapped into the complex via the totally antisymmetric 
three dimensional Lcvi-Civita tensor regarded as thrce-j symbol, i.e. an evaluated three valent 
intertwiner: 



1 



e(e/ (g) ej ® bk) := —i=eijK 
V6 



111 
I J K 



(84) 



Graphically, the result is therefore an ordinary six-j symbol tetrahedral spin network with an 
additional internal vertex, coloured by the evaluated intertwiner ejjK, whose legs, living in the 
adjoint representation, are attached to the three links of the spin network grasped by the source 
derivations corresponding to the wedges of the derived complex T+, i.e to the segments of T. The 
full contribution of the quantum volume of a given vertex v accordingly results in the following 
sum 



grasp 




(85) 



over all sixteen (plus the permutations) possible graspings on the triple of links of the spin network 
which correspond the triple of wedges whose dual segments classically span the volume of the 
tetrahedron dual to the vertex v, i.e over all admissible triples. 

Here, we have introduced a notational simplification that we will use from now on. The formal 
sum and the coefficient K above are entirely defined by the spin network at their right. The 

E 



grasp 



symbol means that we arc summing over all admissible three-graspings of the type defined 
by the spin network. The first term in the sum is the one pictured. For each such term, the 
coefficient X is a function of the spins associated to the grasped faces together with the correct 
sign factor and numerical coefficient. 
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Explicitly, for the above case, we have 



and 



grasp w,w',w" 

1 

ID.SVO 



(86) 



(87) 



The computation of the squared volume of a given vertex v now involves the repeated action 
of Vv on the six-j symbol with source insertions {6j}t,(J) ([77)) associated to v before evaluation 
at J = 0. The result, in terms of spin networks, is the usual tctrahcdral six-j spin network with 
now the insertion of a double trivalent grasping. 

This procedure leads us to a quantisation ambiguity; if at least one of the links belonging 
to the triple of the first three-grasping coincides with one of the links of the second, which is 
automatically satisfied, we have to decide wether the action of the second power of the volume 
grasps "above" or "underneath" the grasping of the first volume action. The two configurations 
yield a different spin network and hence a different number for the quantum measurement of the 
associated tetrahedron's squared volume. To circumvent this difficulty, we propose a natural pre- 
scription consisting in summing over both contributions with equal weight factor, i.e we prescribe 
a symmetrized combination of the two possibilities. Concretely, this means that when we calcu- 
late the squared volume, we first sum over all possible admissible three-graspings (graspings of 
the triple of wedges whose dual segments do not belong to the same triangle). As above, there are 
16 X 3! such terms. For each fixed contribution, we then sum over all the admissible 16 (without 
counting the permutations) second three-grasping possibilities. Out of these terms, there is one 
where the two three-graspings occur on the same triple of links, six where two out of the three 
links are grasped twice and nine contributions where only one link is grasped two times. Our 
prescription is then such that when a term with n links grasped twice occurs in the summation, 
a sum over the 2" possibilities weighted by a factor 1/2" is implemented. Everytime that such 
an ambiguity occurs, the formal notational symbol "^gj-g^^p will contain the implementation of the 
symmetrization procedure. As a result. 
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for each vertex where no fermion is present. Here, 



E --^svA EE' 



(89) 



gras 



and 



^ ■ ^v{jwn jw[: jw" 1 jiU2T jw'^T jw'2) 

2 



(90) 



n ^9niw„ oe(j^je(j^j)e(j^^'). 



where the symmetrization is implied by the symbol SymC^) in the second sum. 
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Fermion loops. We conclude this section by describing the situation along our generic fermionic 
path 7fe . Along his travel round the loop 7^ , the fermion sees a volume contribution Vy living on 
each vertex that he goes trough together with an (internal) area one-form Aje measuring the area 
of each triangle he is travelling through. The component / of the area one-form measured in the 
frame associated to the vertex 1; is calculated by acting with the differential operator 



At, 



1 

12 



E 



^IJK 



5 J J,,, 5J, 



Kw' 



-sgn{w,w'), 



(91) 



on the fermionic seven-j symbol with source insertions {7j(7)}^(J) (|79p living on v. Graphically, 
the result on a given pair of links corresponding to a pair of wedges is the creation of two open 
lines coloured by indices J and K in the adjoint representation. These two lines, together with 
the open line coloured by the index / coming from {7j(7)},f,( J), finally contract trough the ejjK 
intertwining. This is the explicit picture of the soldering of the spin space to the quantum geometry 
of spacetime. 

To complete the picture we also need a volume contribution on the vertex v. Here again the 
procedure implies a double grasping of the fermionic seven-j assigned to the vertices of the fermion 
loop. The result for a given vertex u is the following 



VuAi uvVG 




,7=0 




(92) 



Here, the formal notation means 



and 



K 



E ■■-sy^- E E 



{,jwi 1 jw[ 7 jw2 1 1 3w2 ) 
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(94) 



sgn{wi,w[)sgn{w2,w'^,W2)Q{j^ui)Q{jw^,)Q{jw2)'d{jw2,)'d{jw2,,)- 



The sum is first over the three possible two-graspings of the links corresponding to the wedges 
dual to the segments building up the triangle through which the fermion travels. For each such 
terms there is then a sum over the admissible three-graspings. The Sym symbol ensures that we 
are symmetrizing the appearing ambiguities. 

We are now ready to give a precise definition of the model since all quantities appearing in 
(|57p have been defined. 



3.3.3 The model 

Let us summarise the results and discuss the properties of the model. The proposal made in this 
paper for the description of massive fermionic fields coupled to three dimensional quantum gravity 
consists in an order by order calculation of the partition function (j27p of the theory through an 
inverse mass expansion (|54p. 

The weight of the decorated loop union 7/; = 7p U 7^, p + (7 = fc, of order k € [0, 2n] appearing 
in the a expansion, is given by the following general expression: 



xif) 



^ if 



n E n 



n ^"(^)- 



(95) 
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The no-fermion vertex amplitude B,j is non trivial if and only if the loop ensemble 7^ is built out 
of the union of two loop configurations 7p and 7^ which are such that p and q are smaller than 
n (if one of the two loop configurations saturates all vertices, there is obviously no no-fcrmion 
situation). In this case, the amplitude is defined to be the square volume ([SS]) : Bv{j) ~ . 

The vertex amplitude Ay (7) sitting on the fermion path depends on the type of loop configu- 
ration, decoration considered and also generically varies from one vertex to another. 

For the particular loop configuration 7^. = 7^. U 70 that we have used to build the model, the 
fermionic vertex amplitude has been defined. How does this amplitude generalise to the other 
contributions ? For instance, how to build amplitudes for other types of configurations and other 
types of decorations ? 

We have seen how all the other terms in the sum over loop decorations ([ST]) . i.e all possible 
decorated loops, generated for a fixed set of k vertices could be obtained and we have defined a 
classification of inequivalent loop decorations. Let us pick one general decorated fc-loop G ^^k 
(more precisely, a chosen representative of an equivalence class in Ck/ ^) whose Dirac trace is 
built out of N sequences of the form (jU<jU...(jUW crU^ (J...W a out of our previously established 
classification. Call k ~ m (resp m) the total number of type a : all (rcsp type b: U^cr) like 
patterns appearing in 7^?. Our generic loop hence corresponds to a decorated loop of the form 
71 with TO = 0. Let V-y„ = V jn denote the k vertices traversed by the fermionic loop and 
Va (resp Vb) the subset of vertices associated to the portions of the loop corresponding to the 
type a (resp type b) sub-sequences which are not a starting point for a type b (resp type a) sub- 
sequence. We call Vab and Vba the subset of vertices of V^^ that respectively correspond to an 
end point of a sub-sequence of type a and to an end point of sub-sequence of type b. Hence, 

^VaUVabUVbUVba- 

Taking products of such decorated loops 7, we obtain loop configurations 7. The fermionic 
determinant generates quadratic loop configurations of the form 7^ = 7p U 7g. The value of p and 
q determines three classes of loop configurations. 

Wc now show how the weights associated to all the possible contributions appearing in the 
fermionic loop expansion (j43p can be calculated from the forgoing construction by separately 
considering the three different types of quadratic loop configurations. 

The case p,q < n. If both p and q are strictly smaller than the number of vertices n, we 
generically observe non degenerate vertices. However, the two loop configurations 7p and 7g can 
be supported, partially or totally, by a common set of vertices. We first suppose that 7p H 7^ = 
and consider one component jn C 7p of order r < p. 

Wc now follow a fermion along his travel around the loop starting from a portion of the 
loop associated to a type a sub-sequence belonging to a given previously defined sequence n out 
of the N. Along all the sub-sequence of type a, all vertex amplitudes are given by the following 
vertex functions. Forall v € Va, 



for fc 7^ 2 and 




for the second order case. 
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The fermion now observes a junction between a sub-sequence of type a and a sub-sequence of 
type b yielding a pattern in the trace of the form UU''' . Accordingly, if k ^ 2, all vertices v £ Vab 
carry the following vertex amplitude 



while the k ~ 2 case yields 




because the two group elements U and J7'^ correspond to the same edge: U' = U. The fac- 
tor two comes from the trace formula Tri:,{a^a'^ ) = . Note that the formal sum over 
the graspings here contains a symmetrizing prescription for triple graspings of the same link: 
Egrasp 'S'2/"i(E„,3,to.^,».^' Sym{Y,^,^^^,)) T.w^,w[■ I" ^^is particular case, the sum over the inter- 
twiners in (|95p disappears. It is interesting to remark that, omitting the volume contribution, we 
obtain the expectation value of the area of the triangle dual to the edge on which the fermion is 
travelling. 

Following his trip, the fermion is now going through a type b sub-sequence. All amplitudes for 
k ^ 2 and for vertices v belonging to Vt will be of the form 



^.(7^) = (-1)" E ^ 

grasp 

Finally, the fermions reaches the boundary of a type b sub-sequence to end the sequence n. If 
n ^ N, a, new type a sub-sequence belonging to a new n + I sequence starts. Here two, and only 
two, Pauli matrices are brought together inside the Dirac trace. We accordingly use the relation 
cr^o"^ = ie^j^a^ + rj^'', creating two types of terms associated to the disappearance of the adjacent 
sigma matrices. 

Using the tools developed in this paper, it is not hard to realise that Vw e Vta, the weights are 
given by 



grasp 

Together with the no-fermion vertex amplitude, the above vertex functions define the model 
completely for non-degenerate loop configuration with p,q < n. 

Now, if 7p ^ 0, we observe degenerate vertices traversed by two fermionic lines. If u G Va 
is a degenerate vertex, we obtain the following amplitude forall k ^ 4, and for a given relative 
orientation of the two fermionic lines 




(100) 




(101) 
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or 




(102) 



(103) 



depending on the number of fermionic lines living on the same edge converging towards the 
degenerate vertex. If fc = 4, the two fermionic lines simply close, as in All these amplitudes 
trivially generalise to vertices v £ Vab, Vh, Vba by using the results obtained in the non degenerate 
case. 

The p,q < n case. li p < n and q = n, we necessarily obtain degenerate quadratic loop 
configurations. Some vertices, namely the vertices in the complementary on the 7p configuration, 
however remain non degenerate. If p = q ~ all vertices are degenerate. In all these cases, 
the results developed above generalise yielding the fermionic vertex amplitudes Ay{'-f). Recalling 
that the no-fermion By{'-f) amplitudes are trivial, we have now written down all possible vertex 
amplitudes and the model is consequently defined for all type of loops emerging from the hopping 
parameter expansion. The sum of all these contributions yields the partition function of the 
coupled system. This concludes our presentation of our proposed model of massive fermions 
coupled to three-dimensional quantum gravity. 



4 Conclusion 

In this paper, we have proposed a model describing the coupling of massive fermions to 3d quantum 
gravity by defining the vacuum to vacuum transition amplitude of the theory. We have first defined 
a consistent discretization procedure yielding a concrete meaning to the functional measures. 
Then, we have implemented the Berezin integral over the fermionic fields obtaining the expectation 
value of a functional determinant. This determinant has then been calculated by using a hopping 
parameter expansion yielding a sum over fermionic loops. This finite expansion has enabled us 
to compute the path integral of the theory as we have shown how to calculate the path integral 
order by order in the inverse mass expansion parameter. The integration over simplicial co- frames 
is realised by source derivation on a generalised generating functional where the integration over 
discretized connections is implemented. It is interesting to note that a two-dimensional slicing of 
one of the contributions to the above spinfoam model associated to a given loop, non co-planar 
with the slicing, yields a spin network with an open end coloured by the spinor representation, in 
accordance with loop quantum gravity. The model is rich enough to take into account the effects 
of quantum torsion and of non trivial curvature and highlights the subtle relation between spinors 
and quantum geometry. 

Many issues remain open. First, the emergence of quantisation ambiguities due to the presence 
of high order polynomials in the simplicial triad field in the theory is clearly a problem. We have 
prescribed a coherent procedure but one would need to show that the result does not depend on 
the prescription used to complete the picture. 

Second, unlike the point particle case, fields break the topological invariance of three dimen- 
sional gravity. The exactitude of the model is hence questionable, even if the spinfoam explicitly 
derives from a classical action. This leads to two questions. First, do the asymptotics {j oo) 
of the vertex amplitudes reproduce (a simple function of) the classical Rcgge action coupled to 
fermions ? Second, and most importantly : what is the behaviour of the model under refinement 
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of the triangulation? Since the finiteness of our loop expansion is due to the cut-off on the number 
of degrees of freedom induced by the triangulation, this issue is of primary interest, especially be- 
cause of the large value of the expansion parameter. An alternative to these issues is to consider 
the GFT approach |21| to implement a sum over topologies and triangulations needed when one 
works with non-topological theories. 

Third, the model obviously diverges in the infrared as for pure (non properly gaugcd-fixed 
|18j ) gravity because of the infinite sum over representations. A possible regularization scheme 
to explore consists in incorporating a positive cosmological constant A in the theory leading to 
a modified Turaev-Viro model [U] based on the representation theory of the real compact Hopf 
algebra ilg(su(2)) with the deformation parameter q at root of unity related to the cosmological 
constant A via q ~ e*'p^. Since there are only a finite number of representations of ilq(sM(2)) 
(with non-vanishing g-dimension) at root of unity, this procedure would cut of the infinite sum 
over representations and yield an infrared finite model. 

Finally, one would need to analyse the well known fcrmion doubling problem appearing in 
lattice-like theories incorporating fermions. Does the spinfoam generate doublers ? Could we 
modify the action to suppress the extra species ? 

It is interesting to underline the crucial differences between the coupling of fields and point 
particle to quantum gravity. First, the effects of fields on geometry are non-localised, in the 
sense that matter fields create non trivial vertex amplitudes (squared volumes for instance) on all 
vertices of the dual complex. These effects are more difficult to control and more complicated than 
in the point particle case where curvature and torsion are non trivial only around the particle's 
worldline (note however that this could be due to the mass term). Second, as discussed above, 
fields break the topological invariance of three dimensional gravity for the same reasons, whereas 
particles do not. 

To conclude, we discuss the perspectives. The first natural extention seems to be the gen- 
eralisation to four spacetime dimensions. At first, one might expect the coupling of fermions to 
the Barrett-Crane model of four dimensional quantum gravity to be problematic. As remarked in 
pS] , the reason is that the kinetic term of the standard Dirac action in a d-dimensional curved 
spacetime couples to the d~\ exterior power of the d-bein. Hence, in even dimensions, fermions 
are coupled to Palatini gravity via an odd power of the gravitational field not convertible into 
any power of the d — 2-form B field of BF theory. However, if one uses the Capovilla, Dell, 
Jacobson, Mason chiral action [46] of fermions coupled to two-form gravity, the fermions indeed 
couple to the B field and one could imagine such an action as an interesting starting point for a 
four dimensional construction 0. Next, let us comment on the Lorentzian framework. Because 
of the non-compacity of the associated symmetry group, the unitary, irreducible representations, 
appearing in the character expansion of the delta function forcing the discretized curvature to be 
trivial around the faces of the spinfoam, would be infinite dimensional. Since the fermions live in 
finite dimensional representations, we may expect a problem in the integration over the discretized 
connections since the space of non- vanishing intcrtwincrs between the finite and the unitary rep- 
resentations is empty. Hence, a generalisation of the present work to Lorentzian signatures seems 
a priori non-trivial. Finally, inspired by the recent series of papers devoted to the computation 
of the graviton propagator [47| , it would be interesting to compute the two point function of the 
model. The techniques developed in this paper should be easily generalisable to the construction 
of the fermion propagator in three-dimensional spinfoam quantum gravity. 
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A Appendix : Geometrical and algebraic setting 



This appendix is devoted to review the geometrical and algebraic structures underlying the cou- 
pling of fcrmions to a three dimensional dynamical spacetime. We first define the Clifford algebra 
of interest in this paper, its associated spin group and spinors before introducing the notion of 
spinor bundles and spin connections. Most of the results reviewed here are inspired by references 



A.l Fermions in flat spacetime 

Flat Euclidean spacetime spinors arc related to, 

A. 1.1 The Clifford algebra C(3) and the spin group Spin{3) 

The Clifford algebra C(3,0) = C(3) is the real associative algebra generated by a unit I and the 
three symbols 7^,7^,7'^ satisfying 



7V+7V = 2^", where ^^•^=1^ ^I^j^^'^'^ (104) 

A possible vectorial basis of C(3) is provided by {1, 7^, 7^7"^, J^J'^'y^}i<j<K- The Clifford algebra 
C(3) is therefore a linear space of dimension eight isomorphic to the exterior algebra A(R*^) as a 
vector space and to the space of two by two complex matrices M2 (C) as a real associative algebra. 
The later isomorphism provides us with a faithful representation of C(3) on C^: 

p:C{3) EndiC^) (105) 

7' ^ P(7')=^', 



in terms of the Pauli matrices cr^, I = 1,2, 3, 



1 \ 2 r 2. f 3 , ..^ / 1 



known as the Pauli algebra. The unit is mapped onto the unit two dimensional matrix. Let 
Co (3) be the even subalgebra of C(3), i.e the subalgebra linearly generated by products of an 
even number of generators 7 (i.e two generators) and let Ci(3) denote the vector subspace of C(3) 
linearly generated by the elements 

The group Spm{3) can be explicitly realized as a subset of the Clifford algebra C(3). To see 
this, let V denote a real vector space of dimension three endowed with a non-degenerated bilinear 
form g and denote by {e/}/ a g-orthonormal basis of V: g{ei, ej) = rjij. It is possible to embed 
the algebraic dual V* = Span{e^}/ (endowed with the metric with components 77^'^) of V in C(3) 
through the linear map 



j:V* Ci(3)cC(3) (106) 

V = vie^ — > 7(u) — vi^^ 1^. 

Consider the subset of elements s e C(3) that are invertible and such that 

yveV*, xis)hiv)) = sj{v)s-' eCii'S). (107) 



32 



The set of such elements s forms a group called the Clifford group T. Note that the map x(s) : 
Ci(3) Ci(3) provides us with a linear transformation of V because V, V* and jiV*) = Ci(3) 
are isomorphic as vector spaces. The first isomorphism is due to the fact that the bilinear form g 
is non-degenerate and the second one is given by (|106p . More precisely, x(^) is a-n isometry of the 
bilinear form g, i.e an orthogonal transformation belonging to 0(3). The representation of the 
Clifford group given by x : T ^ GL{V*) is not faithful because, for instance, s and as, a G M*, yield 
the same rotation through x- A standart way of distinguishing between s and as is to introduce 
a normalisation. Spin{3) is then derived from the Clifford group by using this normalisation 
condition. Denote by a "bar" the involution of C(3) defined by 7^72. ..7P = 7^'...7^7-'^. One can 
now define the spin group as the following subgroup of F 

Spin{3) ^{sern Co(3) such that | ss \= ss = 1}, (108) 

where | . | denotes the absolute value. One can show that x{Spin{3)) ~ 5*0(3) and that the map 
X '■ Spin{3) — + 5*0(3) is a surjective two to one group morphism which implies that G = Spin{3) ~ 
SU{2) is the universal double cover of G = 5*0(3) because Spin{3) is simply connected. We can 
equivalently obtain Spin{3) in its fundamental or defining representation as a two by two complex 
matrix group through the isomorphism p. Indeed, in M2(C), the above definition is equivalent 
to defining the group Spin{3) as the group of invertible 2x2 matrices p(s) = A, such that there 
exists a 50(3) linear transformation A{A) : V* V* , satisfying: 

Aaie^)A-^ ^ a{A^j{A)e-^), and det A = 1, (109) 

where a = po ^ : V* ^ M2{C) and the constraint on the determinant forces the element ,s G F to 
belong to the pair subalgebra Co (3) and to satisfy the normalisation condition. We here recognise 
the fact that the adjoint action of SU{2) on its Lie algebra su(2) ~ is a rotation of R'^. 
Physically, we can read this equation as stating that an orthogonal transformation in spacetime 
V must always be accompanied by a change of basis in spin space V = for the spin rotation 
a{v),v € V, to be independent of the frame in which the vector v is expressed. This equation 
also contains information about the profound difference between spinors and vectors, namely that 
spinors are not invariant under transformations corresponding to non trivial elements of the center 
of SU (2). Indeed, note that a 27r radian rotation changes the sign of the spin matrix A of equation 
(|109|) keeping the associated spacetime rotation unchanged. Therefore a spinor on which the spin 
matrix acts will feel the difference between a 2n and a At: rotation whereas a vector will not. 

A. 1.2 Spinors 

We can now define a fermion. Essentially, a fermion ip is a physical field whose local excitations 
yield half integer spin particles and whose dynamics are governed by the Dirac Lagrangian. This 
Lagrangian is built out of complex linear combinations of the fields. Hence, the mathematical 
definition of a classical spin one-half fermion relies on complexified Clifford algebras. 

A Dirac spinor is an element of a complex vector space V on which the complexified even 
Clifford subalgebra Co (3)"' = Co (3) (E)r C is irreducibly represented. Because of the isomorphisms 
of complex algebras Co(3)''' ~ H*' ~ A/2(C), where HI denotes the quaternions, we see that this 
vector space, also called Clifford module, is given by V = C^. Weyl spinors do not exist here 
because Co(3)^ is simple: the algebra cannot be broken into a left and a right part by the use 
of the chirality operator = icr^a^a^. Moreover, (true) Majorana spinors do not exist either 
because the representation V is not of real type. In fact, the representation V is quaternionic. This 
means that V can be endowed with a quaternionic structure, i.e a C-linear[3 function O : V ^ V, 
mapping V onto its complex conjugate vector space V, satisfying C oC ~ — idy, where O : V ^ V 

^^We could also define a quaternionic structure as an anti-linear, anti-involutivc map C : V — ► V. The two 
definitions are equivalent. 



33 



denotes the complex conjugate map to C [49] . This structure aUows us to define a modified version 
of Majorana spinors called symplectic Majorana fermions. 

Because we have the inclusion (|108|) Spin{3) C Co (3), the Dirac spinors -0 e V also yield a 
representation of Spin{3). In fact, this representation V is the fundamental j — 1/2 representation 
of Spin{3), where j G N/2 denotes the spin labelling the irreducible, unitary representations of the 

compact group. Let us introduce some notation used troughout the paper. If V= C^^+^ = ^{6^ \ 

1/2 

a = —j, j}, denotes the left S'pm(3)-modulc associated to the representation j, we have V := V. 

J j j 

We will note V* the contragredient representation space and {e"}a the basis dual to {ea}a- We 

j j j - j 

will call TT G Aut{Y) the representation matrix acting on V. The representation j of G on V induces 

j j j 

a representation of g on V. We will note tt* g Endiy) the associated representation matrix. The 

representations of G are obtained from those of G. Namely, they are the representations of G that 
do not see the center, i.e since (±I) = (±1)^-', we have Wg € G,^ (A(g)) =Tr (g) if and only if 

In order to define a Lagrangian density involving fermions, we need to be able to construct 
scalar objects with spinors. This implies the definition of a S'pm(3)-invariant scalar product on 
V. In pseudo-ricmannian spacetime signatures, this bilinear form involves the introduction of the 
Dirac conjugate (V', 0) = ip4> = ip^j'^cj) G C, where ^ denotes Hermitian conjugation on C and 7° 
is the gamma matrix (representing the 7 symbol) corresponding to the timelike direction. It is 
necessary to introduce the 7*^ matrix to compensate the transformation property of the gamma 
matrices under hermitian conjugation 7^ t — jO^fi^o order to construct a real action. However, 
in the case at hand the gamma matrices arc all hermitian and it is sufficient to define the Dirac 
conjugate by using the 5pm(3)-invariant hermitian bilinear form on V: </>) = ip^P — ''P^4' € C. 
Equivalently, one could consider a bilinear form pairing V with its dual complex conjugate vector 
space V . 

A. 2 Fermions and spacetime geometry 

We now set up the geometrical structure of the spacetime on which the fermions are evolving. 
A. 2.1 Einstein-Cartan geometry 

Let the three-dimensional differential manifold A/ model our spacetime. Let FM denote the frame 
bundle over M, that is, the principal bundle with base M and structure group Gi(3,R) whose 
typical fiber over each point p G AI consists in the set of all possible basis of the tangent space 
TpM . Consider a g[(3, ]R)-valued linear connection on FM. Our spacetime becomes a linearly 
connected space as tangent vectors over different points of the manifold can be compared. One 
can furthermore endow M with an Euclidean metric i.e a non-degenerate bilinear form on the 
tangent space TpM over each point p of the manifold with purely positive signature, and obtain 
a metric affine space where the metric and the linear connection are independent geometrical 
objects. Consider now the reduction of FAI to the G = S'0(3)-bundle of orthonormal frames V 
characterised by the metric g. In the generic case, the linear connection does not know about this 
reduction as it is independent of the metric. Locally, for a given section (a frame, i.e an origin 
and a set of three linearly independent vectors) e : M V\p {e(p)/}/, / = 1,2,3, or for a 
given coordinate basis {-§-}^^ A* = li 2,3, of TpM, the metric is given by 

g = rjije^ (g) e/ ^ g^^dxf^ dx" . (110) 

If we now suppose that the chosen linear connection is compatible with the bundle reduction (i.e 
that it is reducible) meaning that the restriction of the linear connection to the bundle V is a 
principal connection on V, the connection becomes related to the spacetime metric g. In that case, 
the connection obeys the so-called metricity condition stating that the covariant derivative of the 
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metric with respect to the connection vanishes. It is equivalent to simply considering a so(3)- 
valued connection oj on V. The pull-back by a local section of a such connection is a one- form 
on M taking value in the Lie algebra so(3) of 5*0(3), i.c it is matrix valued. If {Taja, a = 1, 2, 3, 
denotes a vectorial basis of so(3), then in a local chart {U C AI,x^ : U M.^) in the domain of 
the local section, 

u; = {u;'j) = ^{u;-'j) = u:l jdx^ ^u^® Tjjdx^. (Ill) 

A such connection is called metric and ensures that the lengths and angles are conserved by parallel 
transport on M. The data of a manifold M endowed with a metric and a metric (compatible) 
connection is called a Riemann-Cartan structure. We call D the covariant derivative associated 
to the metric connection. 

The curvature R ~ Duo of the connection uj is given by the matrix valued two-form 

R = {R'j) = Du'j = dLo'j + A w5 = \R,ljdx^^ A dx^ (112) 

Next, consider the image in P by a local section of a given point p E M. The result is a basis of 
V = M^, i.e a set of three linearly independent vectors e/, orthonormal for the metric g. Let 
denote the associated dual co-frame. It decomposes into a coordinate basis of the cotangent space: 

= ejj^dx^ . The matrix is an isomorphism of vector space between the tangent space over p 
and V regarded as the vector representation space of SO{3); G TpM,e^{v) = € V, i.e it 
expresses any tangent vector into a locally inertial frame where geodesic trajectories are mapped 
onto straight lines: the triad is the gravitational field. We can now define the soldering one 
form with respect to a given choice of a frame and a coframe e : TM — > TM = V X5o(3) ^ as 
e = e-^ (g) £/ = e^dx^ (® e/. 

The torsion T = De of the connection lo is given by the vector valued two-form 

T' = De^ = de' + uj'ii A e^' = ^T'^^dxf" A dx"" . (113) 

In the same way that we have reduced the metric affine structure down to a Ricmann-Cartan 
structure by imposing the metricity condition, we can recover the standard riemannian framework 
of Einstein's general relativity by restricting the set of metric connections to the torsion free or 
riemannian connections. In fact, there is a unique metric compatible, torsion free linear connection. 
It is called the Levi-Civita connection. The Riemann-Cartan geometry |50| that we are considering 
here is therefore a generalisation of ordinary GR. The idea of extending GR as to include non-zero 
torsion solutions was proposed by Elie Cartan in the early twenties. He suggested, before the 
introduction of the modern concept of spin, to relate spacetimc torsion to an intrinsic angular 
momentum of matter. This theory can be obtained as a gauge theory of the Poincarc group where 
torsion and curvature are respectively regarded as the translational and rotational field strengths. 
As we have seen, E. Cartan was right: fermions indeed couple to spacetimc torsion. 

A. 2. 2 Spin structure, Spinor field 

We now wish to introduce fermions on our Ricmann-Cartan spacetimc. Accordingly , we equip 
the manifold M with an extra structure called a spin structure. A spin structure on (M^g), if it 
existsEl, is a principal G-bundle V with base M and structure group G = Spin{3), such that there 
exists a two-to-onc bundle homomorphism x from V onto the bundle V of orthonormal frames 
Troughout this paper, we have assumed that a spin structure had been chosen. Note that all 

^■^The existence of a spin structure over a manifold AI is related to the vanishing of a given characteristic class 
called the Stiefel-Whitney class. 

^^A bundle homomorphism is a map between bundles preserving the fiber structure, i.e commuting with the 
right action of the structure group on the fiber and compatible with the associated group morphism: here, x{p.g) = 
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principal bundles with compact three dimensional base manifold and compact, connected, simply 
connected structural group are necessarily trivial. Our spin bundle is therefore trivial and we can 
consequently choose a global trivialising section once and forall. 

The spin bundle over M being given, we can define a spinor field living on our spacctime. 
A tangent vector to M can be viewed as an element of the vector bundle associated to the 
orthonormal frame bundle V via the fundamental representation of 5*0(3) on F = R'^. In the 
same way, a spinor field ip is defined as a section of the vector bundle SM associated to the 
principal bundle P trough the representation of Spin{3) on V: SAI = Vx spin{3) V. 

To be able to compare the value of the fermionic fields at different points of spacetime, we 
need a covariant derivative acting on spinors. 

A. 2. 3 Spin connection 

A spin connection is a principal connection on P such that its image under the homomorphism 
X is a metric connection on V. The explicit form of the spin connection is obtained from the 
metric connection through the Lie algebra isomorphism x* between spin(3) and 5o(3). By using 
the defining relation (|109p of Spin{3) at the level of the Lie algebras (i.e by differentiation), we 
obtain a unique correspondence between the generators of spin(3) and those of so(3). In a chart 
with local coordinates x'^, the components of the spin(3)-valued connection one- form uj are given 
in terms of the components of the metric connection : 

<^^i^ -^'^"mj, (114) 

where indices are raised and lowered with the metric i]. In the representation of spin(3) 
induced by the representation p of the Clifford algebra C(3) in terms of the Pauli algebra, the spin 
connection is given by p*{lo^) = —■juij/ciaj 

We can now define the covariant derivative one form Vi/' = V p.ipdx^ of a spinor field G 
T(SM) in the spin connection defined above 

= 5pV7 + (115) 



We also define a covariant derivative for the Dirac conjugate, or co-spinor, ip: 

Vfj,ip = df_tip -ipp^iufj,). (116) 
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